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In an earlier paper, W. Kohn had qualitatively introduced the concept of "nearsightedness" of 
electrons in many-atom systems. It can be viewed as underlying such important ideas as Pauling's 
"chemical bond," "transferability" and Yang's computational principle of "divide and conquer." 
It describes the fact that, for fixed chemical potential, local electronic properties, like the density 
n(r), depend significantly on the effective external potential only at nearby points. Changes of 
that potential, no matter how large, beyond a distance R have limited effects on local electronic 
properties, which rapidly tend to zero as function of R. In the present paper, the concept is first 
sharpened for representative models of uncharged fermions moving in external potentials, followed 
by a discussion of the effects of electron-electron interactions and of perturbing external charges. 



PACS numbers: 71.10.- 



71.15.- 



It is a conventional qualitative wisdom among physi- 
cists and chemists that, in the absence of long range ionic 
interactions, if an atom A in a solid is exchanged for an- 
other atom B, the change of the total energy of the sys- 
tem is largely determined by atoms A and B and their 
near neighbors. Similarly for a molecule M adsorbed on 
a surface. These are qualitative examples of "nearsight- 
edness." 

Understanding the physics and chemistry of large 
molecules and solids would have been practically impos- 
sible if not for the principle of transferability 0, ■ ^ 
is generally accepted that, in the absence of long range 
ionic interactions, large molecules or materials systems 
can be studied and understood one neighborhood at a 
time, without the necessity of studying the entire system 
at once. The useful computational method of "divide and 
conquer" takes advantage of this fact [3j . 

We can argue that Pauling's concept Q of the chemi- 
cal bond has a well defined meaning because, to a good 
approximation, its properties depend only on the relative 
positions of the bonded atoms and their near neighbors 
|5j. Anything beyond them, has little influence on the 
properties of the chemical bond. 

These important concepts, based on decades of em- 
pirical and computational work, point to a property of 
matter which we call "nearsightedness of electronic mat- 
ter (NEM)." It was first realized and introduced by one 
of the authors (WK) in 1996 In retrospect, one can 
find precursors of NEM implicit in many other contexts: 
For example, in the work of Thomas and Fermi in the 
1920's |7j; in the proposal of the Local Density Approxi- 
mation in 1965 || and in Yaniv and Kohn'spaper of 1979, 
"Locality Principle in Wave Mechanics" |9j. It was also 
noticed in Lang and Kohn's extensive work on surfaces 
[Io| and as an element of the concept of "edge electrons" 
in Ref. Q. 

NEM deals with the following scenario: We consider an 
unperturbed system of very many charged or uncharged 
electrons in equilibrium in an external, static potential 



v(r), with chemical potential /x, at T = + . We are 
interested in the effect (for fixed fj.) of a perturbing po- 
tential (change of the external potential) w(r'), of finite 
support (footprint) 12], on a local property at a point 
ro, like the density n(r), when the support of w(r') is 
outside a sphere of radius R, centered at ro (Fig. 1). The 
NEM principle states that, for a given unperturbed sys- 
tem and a given R, the density changes at r , An(r ), 
due to all admissible u>(r'), have a finite maximum mag- 
nitude, An, which, of course, depends on ro, R, and on 
the unperturbed system. From this definition, one can 
see that An(ro, R) decays monotonically as a function of 
R. In this paper we prove, for broad classes of systems, 
that in fact, 



R 



lim An(r Q ,R) = 0, 



(1) 



and expect this to be valid very generally. We shall 
show that, for ordered gapless systems, the decay fol- 
lows power-laws, for ordered gapped systems the decay 
is exponential and for disordered, gapped or ungapped 
systems, the decay is also exponential. 

For a given ro and An, we can solve for R from 
An(ro,R) = An and hence define what we call the 
nearsightedness range R(ro,An). The significance of 
R(r , An) is the following: the density changes at ro due 
to any perturbation, of arbitrary shape and amplitude, 
beyond R(ro,An), cannot exceed An. We can say, an- 
thropomorphically that, to within an accuracy An, the 
particle density n (r) cannot "see" any perturbation w(r') 
beyond the distance R(ro,An); hence our word "near- 
sightedness." 

NEM frequently reminds of other well known and well 
understood concepts but, in fact, is different. To avoid 
"deadly sins," let us present a list of what NEM is not: 

1) NEM is not an aspect of linear or higher order non- 
linear response to external perturbations (but does not 
exclude these). 

2) NEM is not screening of charges, which renders long 



FIG. 1: Schematic of Nearsightedness of Electronic Matter: 
v(r) is the unperturbed external potential, w(r') is the per- 
turbing potential outside a sphere of radius R, which is cen- 
tered on the point of interest ro (see text for details) . 



range Coulomb potentials short range (NEM applies also 
to neutral fermions). 

3) NEM does not apply to systems of few electrons or 
to non-interacting bosons below their condensation tem- 
perature (interacting bosons are beyond the scope of this 
paper). 

4) NEM is not limited to electrons at T = 0+ but 
carries over to finite T, including the classical (high T) 
limit. 

5) NEM is not limited to macroscopically homogeneous 
systems. E.g., it applies to a point r on an interface. 

In this article, we communicate the first quantitative 
results on NEM for 1, 2 and 3D non-interacting, peri- 
odic electrons and preliminary results for non-periodic 
and interacting electrons. We shall see that, no mat- 
ter how complicated or strong w(r') is, far away from 
the perturbation, the change of electron density has a 
universal form, which is completely determined by the 
reflection coefficient, in ID, or elements of the scattering 
matrix, in 2 and 3D, evaluated at certain energies. NEM 
follows from the fact that these coefficients cannot exceed 
a certain upper bound. Based on these asymptotic esti- 
mates, we discuss the nearsightedness range and present 
an application to linear scaling electronic algorithms. 

Non-interacting fermions: We emphasise again 
that NEM, as a general principle, does not require in- 
teractions or screening. It is due to the destructive in- 
terference of density (not wave) amplitudes nj associated 
with the occupied single particle eigenstates ipj . 

One dimension: We begin with a model of ID electrons 
in a periodic potential v(x) with inversion symmetry, at 
T = + . The unperturbed Hamiltonian is [H — 2m = 1] 

Ho = -d 2 /dx 2 +v(x), v(x + b) = v(x). (2) 



FIG. 2: The contours Ci and C2 on the Riemann sheet of 
the last occupied band (assumed odd here). Ci surrounds 
the branch point n and the corresponding branch cut (dashed 
line). C2 intersects the real axis at &f and reappears at — kp. 



We first restrict the perturbing potential w(x) to vanish 
for x > 0. The density change is given by 

An(x) = — f [G B (x, x) - G%(x, x)]dE, (3) 
m J c 

where G E and Ge are the unperturbed and perturbed 
Green's functions, respectively, and C is a contour sur- 
rounding the eigenvalues below /i. The integral in Eq. © 
can be mapped into the complex fc-plane, 

An(x) = 2 l R(k)t(; k (x) 2 dk, x > 0, (4) 
Jc 

where C corresponds to C, R(k) is the reflection coeffi- 
cient from right to left and iphix) are the normalized, 
unperturbed Bloch functions 13]. 

a) Asymptotics: For insulators, we can and shall re- 
strict ourselves to Imfc > and to the first Brillouin 
zone. We denote by k the branch point that connects 
the highest occupied and the lowest unoccupied bands 
[14]] . For the case when w(x) generates no bound states 
in the insulating gap, we can choose C in Eq. (gj to be 
the contour C\ in Fig. 2. Since ^k{x) = u k (x)e lkx , with 
Uk(x) periodic of x, ipk{x) decays exponentially with x 
(~ g-imfcz) f or jj-nfc > an j thg asymptotic behavior 
of An(x) comes from the points of C\ in the immediate 
vicinity of k. Using the behavior of the Bloch functions 
near the branch point [lij ]. we find 

/ 2 \ 1/2 

An{x) -> 2R(k) ( — J s K (x) 2 e~ 2 ^, (5) 

where q = Iulk and s K (x) = [\k— K| 1 ' 4 Ufe(a:)]fe^ K is a real, 
periodic or antiperiodic function, s K (x + b) = ±s K (x), de- 
pending on the band index. w(x) enters in this asymp- 
totic form only through R(k). 

The sign of this An(x) is independent of x and is given 
by the sign of R(k). The exponential decay constant q in 



3 



Eq. J5J) equals that of the Wannier function of the highest 
occupied band 0, 0, , or of the density matrix 0] . 
In the limit of a sufficiently small insulating gap G (in 
units of the width W of the highest occupied band), q = 
\^pm*G, where m* is the effective mass at the top of the 
last occupied band [19j. 

One can show that |-R(/«)|max = lj with the maximum 
taken over all functions w(x) which generate no bound 
states in the insulating gap [l3j |. Thus, An(x) cannot 
exceed an upper bound, independent of w (If there are 
bound states, the asymptotic behavior of An(x), Eq. J5J, 
may change, but NEM remains ^3). 

For metals, we integrate Eq. along the contour C 2 
in Fig. 2. For large x, 



An(x) -> -lm[R{k F )^ kF {xf 
x 



(6) 



the slowly decaying Friedel oscillations [2(j. Again, NEM 
follows from the fact that \R(kp)\ cannot exceed 1, for 
any w(x). 

For perturbed Kronig-Penney models |2lj . we found 
that the asymptotic expressions Eqs. JHJ and © set in 
after one or two lattice parameters. 

b) Nearsightedness range R: For a given An, the near- 
sightedness range R(ro,An) at ro was introduced as 
the smallest distance such that any scalar perturbation 
w(r'), lying entirely outside this range, produces a den- 
sity change at r , An(r ), smaller than An. Fully char- 
acterized, R is a function of r , chemical potential fj, and 
An, and a functional of v(r) 



R = R(r >(r')],u,An). 



(7) 



From this definition, it follows that, at every fixed ro, 
dR(r , An)/dAn < 0. 

To calculate R at a point x$ , we need to simultaneously 
consider perturbing potentials wl,r to the left and right 
of xq. In this case, the density change at xq is given by 
the individual contributions Aul,r(xo) ofuiL,R, Eqs. © 
and ©, plus multiple reflections corrections. In the limit 
when the distance from xo to wl.r is large, these correc- 
tions were found exponentially small for insulators and 
comparable to An^(xo) + An^(xo) for metals. 

For insulators, the asymptotic behavior of R in the 
limit An — * 0, as derived from the upper bound of 
Eq. @ , from the above remark and from a cell-averaging 
of An(x), is 



R(An) 



1 



■In — 

2q Ar 



where 



dx. 



(9) 



In the small gap and tight binding limits, n is completely 
determined by the exponential decay constant q, h — > 
Aqy / 2/n and by h — > 4^/ q/irb, respectively. 



For metals, the upper bound of Eq. JfjJ, with inclu- 
sion of multiple reflections and cell-averaging, leads to 
the following asymptotic behavior 



R(An) 1/An. 



(10) 



c) Disorder: NEM is not limited to periodic v(x). As 
a first orientation, we discuss the effect of adding a small 
random potential 



v r (x) = A 5 (a 



(11) 



to the periodic potential v(x) in Eq. J2J. We have calcu- 
lated the averaged density-density correlation function, 
to the lowest order in A, assuming a random distribution 
of impurities with an average density n^ (cf. 22]). Our 
conclusions are as follows. 

For insulators, the random potential changes the ex- 
ponential decay constant q by 



Aq = —riiX 



: 3(4tt) 2 <■» 
26 7 2 



[s- K (x)s K (x)} 2 dx, (12) 



where 7 = [(E k - E K )/(k ~ k) 1 / 2 ]^. Two thirds of Aq 
are due to the narrowing of the gap G, 



AG= ~ niX2 ^G~~ /Vk(^K(*)] 2 dx, 



(13) 



and one third is due to the fluctuations. 

For metals, the random potential introduces an expo- 
nential decay in the averaged density response, with a 
decay constant 



q = ni\ 2 g(e F ) 



,12tt 



4 r b 



|V>fc F (a;)| 4 (i:r, 



(14) 



where g(ep) is the density of states at the Fermi energy. 

Thus, a small random potential increases R for insula- 
tors and decreases it for metals. We conjecture that this 
remains true for disorder, more generally. 

Higher dimensions: We first consider, as an example, 
2D fermions in a periodic potential with square symme- 
try. We present only the cases when the first band, as- 
sumed simple and isolated (2^| , is partially or completely 
filled. We first restrict the perturbing potential w(x,y) 
to vanish for y > and the periodicity along the x direc- 
tion to be preserved (so that k x remains a good quantum 
number). The density change for y — > 00 is given by [24| 

An(f)->2 J S(k,k')il; R (r)^i(r)dk, (15) 

with k! = (k x ,—ky) and S(k,k') the scattering matrix 
element between k and k' . 

a) Asymptotics: For fixed k x , the analytic structure 
relative to k y of the band energy E^ and Bloch function 
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of the first band is completely analogous to that in 
ID Et and ipr have branch points of order 1 and 

3, respectively, at K y = ±7r + iq(k x ), connecting the first 
band with a higher band (we restrict ourselves to Imk y > 
and to the first Brillouin zone). Their behavior near 
these points is the same as in ID, namely E? behaves as 
a square root and ip? diverges as (k y — K y ) -1 / 4 j^. 

If the band is completely filled and w(r) does not gen- 
erate bound states in the insulating gap, for a given 
k x , the integral over k y in Eq. can be taken over 
a contour surrounding the branch point (as in Fig. 2) 
and its asymptotic behavior can be determined as in the 
ID case. Thus, the integral over k y decays exponen- 
tially, as function of y, with a rate 2q(k x ). There will 
be two values, ±fe£, of k x where q(k x ) reaches its lowest 
value go- The asymptotic behavior of An(r), for y — > 
oo, comes from the immediate vicinity of these points. 
Defining j3 = (d 2 q/d 2 k x ) kx = k o and k = (k°, tt + iq ), 
4 = { k x,-n- i( lo), gives 



An(r) 



4tt 

7V 



Re[S(K , Z' )fa (mdr)}e- 2qa \ (16) 



where <fe (f) ee [\k y - K y \ 1/4 e lkyV tp k {r)}^^ is a quasi- 
periodic function in x and y. Again, w(r) enters in this 
asymptotic form only through S(kq,k' ). 

The exponential decay in Eq. l|16f) is twice as fast as 
the exponential decay, in the y direction, of the den- 
sity matrix or of the Wannier function of the first band. 
\S(kq, k' q )\ = 1 for a hard wall and, in general, we expect 
it to be of order 1. 

If the band is partially filled, the asymptotic behavior 
of An(r) is determined by the two points on the Fermi 
surface, denoted by ko = {k x ,k y ) and k' = (k x ,—k y ), 
where the tangent to the Fermi surface is along the k x - 
direction: 



An(r) 



21m* 



S(k ,k' )^(r)r P (^, (17) 



with r] the curvature of the Fermi surface at these points. 
From the unitarity of the scattering matrix, one can im- 
mediately find that \S(ko, k' )\ max — 1, with the maxi- 
mum taken over all w(r). Thus, the asymptotic density 
change cannot exceed an upper bound, no matter how 
large the perturbing potential is. 
The 3D case is analogous. 

b) Nearsightedness range R: For metals, the simplest 
model is jellium enclosed by a spherical hard wall, for 
which one easily finds 



R(An) 



k F /2An = 2.2r s n/An (2D) 
k 2 F /2irAn = 2.5r s n/An {3D), 



(18) 



where n is the density of the uniform gas and r s the 
Wigner-Seitz radius. For metals in periodic potentials, 



the An dependence remains unchanged, but kp is re- 
placed by a k e ff depending on the band structure and 
the filling. 

For an insulator with square symmetry, we calculated 
An(ro) due to enclosing the point ro in four hard walls 
along the symmetry axes, at a distance D from ro- The 
density change near ro is given by the sum of the changes 
due to each individual, infinitely extended wall, Eq. I|I6(I . 
plus multiple reflection corrections, which were found to 
be exponentially negligible in the limit D — > oo. Similarly 
for a 3D insulator with cubic symmetry. From Eq. JT^J, 
its 3D analog |24| and the previous remark, we find that 
the cell-averaged density change at r becomes less than 
a given An for D > D(An), 



I n 
D(An) -> — In—, 
v ; 2q An 



(19) 



where h can be easily calculated from the band structure. 

Finding an analytic expression of the nearsightedness 
range R for a general 2 or 3D insulator is clearly a next 
to impossible task. However, on the basis of the above 
calculations, the proof of exp onential localization of 2 
and 3D Wannier functions |2jj|,|26j and of ID generalized 
Wannier functions 0,^1, as we ll as our ID result, Eqs. 
(0 and JSJ), we expect results of the following form for 2 
and 3D insulators: 



R(An) 



2q e ff An 



(20) 



where g e // is an exponential decay constant of the density 
matrix and h oc q^ ^ (d = 2, 3). 

c) Disorder: We expect the effects of disorder in 2 and 
3D to be qualitatively similar to those in ID. 

Linear scaling: The CPU time for electronic struc- 
ture calculations of a system consisting of many (N a ) 
atoms grows very rapidly with N a , if the calculations are 
performed for the entire system at once. It has been 
pointed out 0, H3 that the dependence on N a can be 
made linear for large N a , by dividing the system into N s 
suitable sub-systems, where N s oc N a . In Ref. [||, NEM 
was identified as the physical basis of linear scaling. Here, 
we quantify this idea. 

The procedure is illustrated in Fig. 3 (a detailed dis- 
cussion (including self-consistency) can be found in Ref. 
[13 )■ The system of volume V is divided into segments 
V n , with overlapping buffer zones B n . The density n(r) 
and contribution to the total energy of each V n are ob- 
tained from calculations including its buffer zone. For 
a required accuracy An, the thickness b of B n is chosen 
so that, when r is on the boundary of V n , An(r) < An, 
where An(r) is the error due to the hard walls around 
B n . An upper bound for b is given by the maximum of 
R(r, An) on the boundary of V n . Finally, the size of each 
V n is chosen to minimize the total CPU time. 

We exemplify this for "cubic" periodic systems in 1, 2 
and 3D, where the V n are all identical "cubes" of edge size 



5 























Vn 

a 


b 
















/ 

Hard wall 





FIG. 3: The system is divided into smaller volumes V n (9 in 
this example), with buffer zones B n (gray color). 



a. The largest density change occurs at the "corners" of 
V n and comes primarily from the nearest walls. For met- 
als, Eqs. JJJJ, (|17|) and their 3D analog, plus the inclusion 
of the multiple reflections, lead to b — > x^fn/An) 2 '^ 1 ', 
with x = 0.31, 0.91 and 1.0 for 1, 2 and 3D, respectively. 
For a 5% accuracy, b = 6.2, 6.7 and 4.5r s , respectively. 
Similarly, for insulators, b —>■ {2q a )~ l ln[n/2An], with h 
defined in Eqs. @ and (JTSJ. 

The total CPU time is given by i = N 3 r, where N s 
(= V/a d ) is the number of segments, and r (oc (26+a)" d ) 
is the CPU time for the electronic structure calculation 
of one segment plus its buffer zone [v = 2 — 3 for DFT 
[2?| , and higher for other methods] . Minimizing the total 
CPU time with respect to a, we obtain the optimal size, 
a = 2b/ (v — 1). With this optimization, and from our 
estimates of b, we obtain the following dependence of the 
total CPU time on the desired accuracy and the total 
number of atoms: 



t oc N a x 



(lnn^An)^ 1 ^ 



(ungapped) 
(gapped) 



(21) 



For metals, b can be greatly reduced by averaging the 
wall-induced Friedel oscillations over two or more values 
of b. 

Interacting fermions: In considering the response of 
charged fermions to distant disturbances it is necessary 
to distinguish between two cases: 

a) Distant perturbing potentials, w{r'), with |rrj— r'\ > 
R. The simplest description of many-body interaction 
effects is the random phase approximation which, in all 
dimensions, leads to a decrease of R in typical metals but 
an increase of R in typical insulators due to a reduction 
in the gap 



b) Distant perturbing charge densities p(r'). In anal- 
ogy with R(ro,An), we define a charge-nearsightedness 
range, R c (r , An) as the smallest distance such that any 
charge perturbation p(r') lying entirely outside this range 
produces a density change at ro, An(ro), smaller than 
An. 

As is well known, the long range Coulomb potential, 
due to perturbing electric charges, is screened out by 
metallic electrons. Preliminary model calculations for 
metallic electrons, in the Thomas-Fermi approximation, 
indicate that they are charge-nearsighted, i.e. have a fi- 
nite R c . However, charged insulating fermions are "clas- 
sically farsighted," in the sense that, at sufficiently large 
distances, the fermions "see" the classical long range to- 
tal potential / p t (r')/\ro — r'\dr' , were p t is the total per- 
turbing charge density, including depolarization. Thus, 
for example, in metals, replacing a neutral atom or ion 
by another atom or ion has always short range electronic 
consequences, while, in an insulator, ions lead to classical 
long range electronic effects. 
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